    R = zeros(N,N);
    for i = 1:N
        for j = 1:N
            yn = Sigma(i*j);
            R(i,j) = yn(1-yn);
        end
    end
    H = data'*R*data;   %(N*D)'*(N*N)*(N*D)